Plant Gender Affects Soil Fungal Microbiota Associated with Welwitschia mirabilis, an Unusual Desert Gymnosperm.

In a recent study, we found a distinct soil bacterial community associated with male and female plants of the desert gymnosperm Welwitschia mirabilis. In this subsequent study, we also found that the soil fungal community associated with Welwitschia differs between male and female plants, and between unvegetated areas and the soil under plants. Site location, pH, and soil moisture also had an important influence on the composition of the fungal community. A number of Ascomycota and Chytrid species were found to be distinct indicators of male and female plants, respectively, but there was no overall difference at the phylum level or in terms of diversity. The unvegetated areas between plants also differed in terms of several Ascomycota OTUs. Network connectivity of the fungal communities was found to be higher under both male and female Welwitschia plants than in unvegetated control areas. As with the bacterial community, it is unclear what processes produce the gender-distinct fungal community, and also the more general plant-associated community, and also what the effects on the biology of the plants are. One possibility behind the gender-related difference in fungal community is that there are differences in the production of pollen or nectar between the two plant genders, affecting the below-ground soil community.

Projection of Gut Microbiome Pre- and Post-Bariatric Surgery To Predict Surgery Outcome.

Bariatric surgery is often the preferred method to resolve obesity and diabetes, with ∼800,000 cases worldwide yearly and high outcome variability. The ability to predict the long-term body mass index (BMI) change following surgery has important implications for individuals and the health care system in general. Given the tight connection between eating habits, sugar consumption, BMI, and the gut microbiome, we tested whether the microbiome before any treatment is associated with different treatment outcomes, as well as other intakes (high-density lipoproteins [HDL], triglycerides, etc.). A projection of the gut microbiome composition of obese (sampled before and after bariatric surgery) and lean patients into principal components was performed, and the relation between this projection and surgery outcome was studied. The projection revealed three different microbiome profiles belonging to lean, obese, and obese individuals who underwent bariatric surgery, with the postsurgery microbiome more different from the lean microbiome than the obese microbiome. The same projection allowed for a prediction of BMI loss following bariatric surgery, using only the presurgery microbiome. The microbial changes following surgery were an increase in the relative abundance of Proteobacteria and Fusobacteria and a decrease in Firmicutes. The gut microbiome can be decomposed into main components depicting the patient's development and predicting in advance the outcome. Those may be translated into the better clinical management of obese individuals planning to undergo metabolic surgery. IMPORTANCE BMI and diabetes can affect the gut microbiome composition. Bariatric surgery has large variabilities in the outcome. The microbiome was previously shown to be a good predictor for multiple diseases. We analyzed here the gut microbiome before and after bariatric surgery and showed the following. (i) The microbiome before surgery can be used to predict surgery outcomes. (ii) The postsurgery microbiome drifts further away from the lean microbiome than the microbiome of the presurgery obese patients. These results can lead to a microbiome-based presurgery decision whether to perform surgery.

Resistance to paclitaxel is associated with a variant of the gene BCL2 in multiple tumor types.

Paclitaxel, the most commonly used form of chemotherapy, is utilized in curative protocols in different types of cancer. The response to treatment differs among patients. Biological interpretation of a mechanism to explain this personalized response is still unavailable. Since paclitaxel is known to target BCL2 and TUBB1, we used pan-cancer genomic data from hundreds of patients to show that a single-nucleotide variant in the BCL2 sequence can predict a patient's response to paclitaxel. Here, we show a connection between this BCL2 genomic variant, its transcript structure, and protein abundance. We demonstrate these findings in silico, in vitro, in formalin-fixed paraffin-embedded (FFPE) tissue, and in patient lymphocytes. We show that tumors with the specific variant are more resistant to paclitaxel. We also show that tumor and normal cells with the variant express higher levels of BCL2 protein, a phenomenon that we validated in an independent cohort of patients. Our results indicate BCL2 sequence variations as determinants of chemotherapy resistance. The knowledge of individual BCL2 genomic sequences prior to the choice of chemotherapy may improve patient survival. The current work also demonstrates the benefit of community-wide, integrative omics data sources combined with in-lab experimentation and validation sets.

Evaluating methods for Avian avulavirus-1 whole genome sequencing.

BACKGROUND: Avian avulavirus-1 (AAvV-1, previously Newcastle Disease Virus) is responsible for poultry and wild birds' disease outbreaks. Numerous whole genome sequencing methods were reported for this virus. These methods included cloning, specific primers amplification, shotgun PCR approaches, Sequence Independent Single Primer Amplification and next generation sequencing platform kits. METHODS: Three methods were used to sequence 173 Israeli Avian avulavirus-1 field isolates and one vaccine strain (VH). The sequencing was performed on Proton and Ion Torrent Personal Genome Machine and to a lesser extent, Illumina MiSeq and NextSeq sequencers. Target specific primers (SP) and Sequence Independent Single Primer Amplification (SISPA) products sequenced via the Ion torrent sequencer had a high error rate and truncated genomes. All the next generation sequencing platform sequencing kits generated high sequence accuracy and near-complete genomic size. RESULTS: A high level of mutations was observed in the intergenic regions between the avian avulavirus-1 genes. Within genes, multiple regions are more mutated than the Fusion region currently used for typing. CONCLUSIONS: Our findings suggest that the whole genome sequencing by the Ion torrent sequencing kit is sufficient. However, when higher fidelity is desired, the Illumina NextSeq and Proton torrent sequencing kits were found to be preferable.

Evaluating methods for Avian avulavirus-1 whole genome sequencing.

BACKGROUND: Avian avulavirus-1 (AAvV-1, previously Newcastle Disease Virus) is responsible for poultry and wild birds' disease outbreaks. Numerous whole genome sequencing methods were reported for this virus. These methods included cloning, specific primers amplification, shotgun PCR approaches, Sequence Independent Single Primer Amplification and next generation sequencing platform kits. METHODS: Three methods were used to sequence 173 Israeli Avian avulavirus-1 field isolates and one vaccine strain (VH). The sequencing was performed on Proton and Ion Torrent Personal Genome Machine and to a lesser extent, Illumina MiSeq and NextSeq sequencers. Target specific primers (SP) and Sequence Independent Single Primer Amplification (SISPA) products sequenced via the Ion torrent sequencer had a high error rate and truncated genomes. All the next generation sequencing platform sequencing kits generated high sequence accuracy and near-complete genomic size. RESULTS: A high level of mutations was observed in the intergenic regions between the avian avulavirus-1 genes. Within genes, multiple regions are more mutated than the Fusion region currently used for typing. CONCLUSIONS: Our findings suggest that the whole genome sequencing by the Ion torrent sequencing kit is sufficient. However, when higher fidelity is desired, the Illumina NextSeq and Proton torrent sequencing kits were found to be preferable.

Draft Genome Sequence of the Suttonellaornithocola Bacterium.

We report here the draft genome sequence of the Suttonella ornithocola bacterium. To date, this bacterium, found in birds, passed only phylogenetic and phenotypic analyses. To our knowledge, this is the first publication of the Suttonella ornithocola genome sequence. The genetic profile provides a basis for further analysis of its infection pathways.

Draft Genome Sequence of a New Bacterium Named Rappaport israeli, gen. nov., sp. nov., Isolated from a Blood Culture.

Here, we report the draft genome sequence of a Gram-negative microbe found in a blood culture (B08008) from a patient. The organism was proposed to be from a new unknown genus and species. This publication will increase worldwide microbial knowledge and may improve microbial identification and antibiotic treatment for patients.

SELMAP - SELEX affinity landscape MAPping of transcription factor binding sites using integrated microfluidics.

Transcription factors (TFs) alter gene expression in response to changes in the environment through sequence-specific interactions with the DNA. These interactions are best portrayed as a landscape of TF binding affinities. Current methods to study sequence-specific binding preferences suffer from limited dynamic range, sequence bias, lack of specificity and limited throughput. We have developed a microfluidic-based device for SELEX Affinity Landscape MAPping (SELMAP) of TF binding, which allows high-throughput measurement of 16 proteins in parallel. We used it to measure the relative affinities of Pho4, AtERF2 and Btd full-length proteins to millions of different DNA binding sites, and detected both high and low-affinity interactions in equilibrium conditions, generating a comprehensive landscape of the relative TF affinities to all possible DNA 6-mers, and even DNA10-mers with increased sequencing depth. Low quantities of both the TFs and DNA oligomers were sufficient for obtaining high-quality results, significantly reducing experimental costs. SELMAP allows in-depth screening of hundreds of TFs, and provides a means for better understanding of the regulatory processes that govern gene expression.

Desert Perennial Shrubs Shape the Microbial-Community Miscellany in Laimosphere and Phyllosphere Space.

Microbial function, composition, and distribution play a fundamental role in ecosystem ecology. The interaction between desert plants and their associated microbes is expected to greatly affect their response to changes in this harsh environment. Using comparative analyses, we studied the impact of three desert shrubs, Atriplex halimus (A), Artemisia herba-alba (AHA), and Hammada scoparia (HS), on soil- and leaf-associated microbial communities. DNA extracted from the leaf surface and soil samples collected beneath the shrubs were used to study associated microbial diversity using a sequencing survey of variable regions of bacterial 16S rRNA and fungal ribosomal internal transcribed spacer (ITS1). We found that the composition of bacterial and fungal orders is plant-type-specific, indicating that each plant type provides a suitable and unique microenvironment. The different adaptive ecophysiological properties of the three plant species and the differential effect on their associated microbial composition point to the role of adaptation in the shaping of microbial diversity. Overall, our findings suggest a link between plant ecophysiological adaptation as a "temporary host" and the biotic-community parameters in extreme xeric environments.

Identification of novel transcriptional regulators of Zat12 using comprehensive yeast one-hybrid screens.

To appropriately acclimate to environmental stresses, plants have to rapidly activate a specific transcriptional program. Yet, the identity and function of many of the transcriptional regulators that mediate early responses to abiotic stress stimuli is still unknown. In this work we employed the promoter of the multi-stress-responsive zinc-finger protein Zat12 in yeast one-hybrid (Y1H) screens to identify early abiotic stress-responsive transcriptional regulators. Analysis of Zat12 promoter fragments fused to luciferase underlined an approximately 200 bp fragment responsive to NaCl and to reactive oxygen species (ROS). Using these segments and others as baits against Y1H control or stress Arabidopsis prey libraries, we identified 15 potential Zat12 transcriptional regulators. Among the prominent proteins identified were known transcription factors including bZIP29 and ANAC91 as well as unknown function proteins such as a homolog of the human USB1, a U6 small nuclear RNA (snRNA) processing protein, and dormancy/auxin-associated family protein 2 (DRM2). Altered expression of Zat12 during high light stress in the knockout mutants further indicated the involvement of these proteins in the regulation of Zat12. Using a state of the art microfluidic approach we showed that AtUSB1 and DRM2 can specifically bind dsDNA and were able to identify the preferred DNA-binding motif of all four proteins. Overall, the proteins identified in this work provide an important start point for charting the earliest signaling network of Zat12 and of other genes required for acclimation to abiotic stresses.

Reduced levels of protein recoding by A-to-I RNA editing in Alzheimer's disease.

Adenosine to inosine (A-to-I) RNA editing, catalyzed by the ADAR enzyme family, acts on dsRNA structures within pre-mRNA molecules. Editing of the coding part of the mRNA may lead to recoding, amino acid substitution in the resulting protein, possibly modifying its biochemical and biophysical properties. Altered RNA editing patterns have been observed in various neurological pathologies. Here, we present a comprehensive study of recoding by RNA editing in Alzheimer's disease (AD), the most common cause of irreversible dementia. We have used a targeted resequencing approach supplemented by a microfluidic-based high-throughput PCR coupled with next-generation sequencing to accurately quantify A-to-I RNA editing levels in a preselected set of target sites, mostly located within the coding sequence of synaptic genes. Overall, editing levels decreased in AD patients' brain tissues, mainly in the hippocampus and to a lesser degree in the temporal and frontal lobes. Differential RNA editing levels were observed in 35 target sites within 22 genes. These results may shed light on a possible association between the neurodegenerative processes typical for AD and deficient RNA editing.

Fmrp Interacts with Adar and Regulates RNA Editing, Synaptic Density and Locomotor Activity in Zebrafish.

Fragile X syndrome (FXS) is the most frequent inherited form of mental retardation. The cause for this X-linked disorder is the silencing of the fragile X mental retardation 1 (fmr1) gene and the absence of the fragile X mental retardation protein (Fmrp). The RNA-binding protein Fmrp represses protein translation, particularly in synapses. In Drosophila, Fmrp interacts with the adenosine deaminase acting on RNA (Adar) enzymes. Adar enzymes convert adenosine to inosine (A-to-I) and modify the sequence of RNA transcripts. Utilizing the fmr1 zebrafish mutant (fmr1-/-), we studied Fmrp-dependent neuronal circuit formation, behavior, and Adar-mediated RNA editing. By combining behavior analyses and live imaging of single axons and synapses, we showed hyperlocomotor activity, as well as increased axonal branching and synaptic density, in fmr1-/- larvae. We identified thousands of clustered RNA editing sites in the zebrafish transcriptome and showed that Fmrp biochemically interacts with the Adar2a protein. The expression levels of the adar genes and Adar2 protein increased in fmr1-/- zebrafish. Microfluidic-based multiplex PCR coupled with deep sequencing showed a mild increase in A-to-I RNA editing levels in evolutionarily conserved neuronal and synaptic Adar-targets in fmr1-/- larvae. These findings suggest that loss of Fmrp results in increased Adar-mediated RNA editing activity on target-specific RNAs, which, in turn, might alter neuronal circuit formation and behavior in FXS.

Synthesis and evaluation of a photoresponsive quencher for fluorescent hybridization probes.

Nowadays, most nucleic acid detections using fluorescent probes rely on quenching of fluorescence by energy transfer from one fluorophore to another or to a non-fluorescent molecule (quencher). The most widely used quencher in fluorescent probes is 4-((4-(dimethylamino)phenyl)azo)benzoic acid (DABCYL). We targeted a nucleoside-DABCYL analogue which could be incorporated anywhere in an oligonucleotide sequence and in any number, and used as a quencher in different hybridization sensitive probes. Specifically, we introduced a 5-(4-((dimethylamino)phenyl)azo)benzene)-2'-deoxy-uridine (dU(DAB)) quencher. The photoisomerization and dU(DAB)'s ability to quench fluorescein emission have been investigated. We incorporated dU(DAB) into a series of oligonucleotide (ON) probes including strand displacement probes, labeled with both fluorescein (FAM) and dU(DAB), and TaqMan probes bearing one or two dU(DAB) and a FAM fluorophore. We used these probes for the detection of a DNA target in real-time PCR (RT-PCR). All probes showed amplification of targeted DNA. A dU(DAB) modified TaqMan RT-PCR probe was more efficient as compared to a DABCYL bearing probe (93% vs. 87%, respectively). Furthermore, dU(DAB) had a stabilizing effect on the duplex, causing an increase in Tm up to 11 °C. In addition we showed the photoisomerisation of the azobenzene moiety of dU(DAB) and the dU(DAB) triply-labeled oligonucleotide upon irradiation. These findings suggest that dU(DAB) modified probes are promising probes for gene quantification in real-time PCR detection and as photoswitchable devices.

Drosophila TRF2 is a preferential core promoter regulator.

Transcription of protein-coding genes is highly dependent on the RNA polymerase II core promoter. Core promoters, generally defined as the regions that direct transcription initiation, consist of functional core promoter motifs (such as the TATA-box, initiator [Inr], and downstream core promoter element [DPE]) that confer specific properties to the core promoter. The known basal transcription factors that support TATA-dependent transcription are insufficient for in vitro transcription of DPE-dependent promoters. In search of a transcription factor that supports DPE-dependent transcription, we used a biochemical complementation approach and identified the Drosophila TBP (TATA-box-binding protein)-related factor 2 (TRF2) as an enriched factor in the fractions that support DPE-dependent transcription. We demonstrate that the short TRF2 isoform preferentially activates DPE-dependent promoters. DNA microarray analysis reveals the enrichment of DPE promoters among short TRF2 up-regulated genes. Using primer extension analysis and reporter assays, we show the importance of the DPE in transcriptional regulation of TRF2 target genes. It was previously shown that, unlike TBP, TRF2 fails to bind DNA containing TATA-boxes. Using microfluidic affinity analysis, we discovered that short TRF2-bound DNA oligos are enriched for Inr and DPE motifs. Taken together, our findings highlight the role of short TRF2 as a preferential core promoter regulator.

The hnRNP F/H homologue of Trypanosoma brucei is differentially expressed in the two life cycle stages of the parasite and regulates splicing and mRNA stability.

Trypanosomes are protozoan parasites that cycle between a mammalian host (bloodstream form) and an insect host, the Tsetse fly (procyclic stage). In trypanosomes, all mRNAs are trans-spliced as part of their maturation. Genome-wide analysis of trans-splicing indicates the existence of alternative trans-splicing, but little is known regarding RNA-binding proteins that participate in such regulation. In this study, we performed functional analysis of the Trypanosoma brucei heterogeneous nuclear ribonucleoproteins (hnRNP) F/H homologue, a protein known to regulate alternative splicing in metazoa. The hnRNP F/H is highly expressed in the bloodstream form of the parasite, but is also functional in the procyclic form. Transcriptome analyses of RNAi-silenced cells were used to deduce the RNA motif recognized by this protein. A purine rich motif, AAGAA, was enriched in both the regulatory regions flanking the 3' splice site and poly (A) sites of the regulated genes. The motif was further validated using mini-genes carrying wild-type and mutated sequences in the 3' and 5' UTRs, demonstrating the role of hnRNP F/H in mRNA stability and splicing. Biochemical studies confirmed the binding of the protein to this proposed site. The differential expression of the protein and its inverse effects on mRNA level in the two lifecycle stages demonstrate the role of hnRNP F/H in developmental regulation.

Functional analysis of noncoding RNAs in trypanosomes: RNA walk, a novel approach to study RNA-RNA interactions between small RNA and its target.

The recent discovery of thousands of small noncoding RNAs (ncRNAs), in many different organisms, has led to the need for methods to study their function. One way to help understand their function is to determine what other RNAs interact with the ncRNAs. We have developed a novel method to investigate the RNA-RNA interactions between a small RNA and its target that we termed "RNA walk." The method is based on UV-induced AMT cross-linking in vivo followed by affinity selection of the hybrid molecules and mapping the intermolecular adducts by RT-PCR. Domains carrying the cross-linked adducts are less efficiently amplified than domains that are not cross-linked. Real-time PCR is used to quantify the results. Further mapping of the interactions is performed by primer extension to determine the exact cross-linked adduct.

A comparative genome-wide study of ncRNAs in trypanosomatids.

BACKGROUND: Recent studies have provided extensive evidence for multitudes of non-coding RNA (ncRNA) transcripts in a wide range of eukaryotic genomes. ncRNAs are emerging as key players in multiple layers of cellular regulation. With the availability of many whole genome sequences, comparative analysis has become a powerful tool to identify ncRNA molecules. In this study, we performed a systematic genome-wide in silico screen to search for novel small ncRNAs in the genome of Trypanosoma brucei using techniques of comparative genomics. RESULTS: In this study, we identified by comparative genomics, and validated by experimental analysis several novel ncRNAs that are conserved across multiple trypanosomatid genomes. When tested on known ncRNAs, our procedure was capable of finding almost half of the known repertoire through homology over six genomes, and about two-thirds of the known sequences were found in at least four genomes. After filtering, 72 conserved unannotated sequences in at least four genomes were found, 29 of which, ranging in size from 30 to 392 nts, were conserved in all six genomes. Fifty of the 72 candidates in the final set were chosen for experimental validation. Eighteen of the 50 (36%) were shown to be expressed, and for 11 of them a distinct expression product was detected, suggesting that they are short ncRNAs. Using functional experimental assays, five of the candidates were shown to be novel H/ACA and C/D snoRNAs; these included three sequences that appear as singletons in the genome, unlike previously identified snoRNA molecules that are found in clusters. The other candidates appear to be novel ncRNA molecules, and their function is, as yet, unknown. CONCLUSIONS: Using comparative genomic techniques, we predicted 72 sequences as ncRNA candidates in T. brucei. The expression of 50 candidates was tested in laboratory experiments. This resulted in the discovery of 11 novel short ncRNAs in procyclic stage T. brucei, which have homologues in the other trypansomatids. A few of these molecules are snoRNAs, but most of them are novel ncRNA molecules. Based on this study, our analysis suggests that the total number of ncRNAs in trypanosomatids is in the range of several hundred.

Establishment of an in vitro trans-splicing system in Trypanosoma brucei that requires endogenous spliced leader RNA.

In trypanosomes a 39 nucleotide exon, the spliced leader (SL) is donated to all mRNAs from a small RNA, the SL RNA, by trans-splicing. Since the discovery of trans-splicing in trypanosomes two decades ago, numerous attempts failed to reconstitute the reaction in vitro. In this study, a crude whole-cell extract utilizing the endogenous SL RNA and synthetic tubulin pre-mRNA were used to reconstitute the trans-splicing reaction. An RNase protection assay was used to detect the trans-spliced product. The reaction was optimized and shown to depend on ATP and intact U2 and U6 snRNPs. Mutations introduced at the polypyrimidine tract and the AG splice site reduced the reaction efficiency. To simplify the assay, RT-PCR and quantitative real-time PCR assays were established. The system was used to examine the structural requirements for SL RNA as a substrate in the reaction. Interestingly, synthetic SL RNA assembled poorly to its cognate particle and was not utilized in the reaction. However, SL RNA synthesized in cells lacking Sm proteins, which is defective in cap-4 modification, was active in the reaction. This study is the first step towards further elucidating the mechanism of trans-splicing, an essential reaction which determines the trypanosome transcriptome.

'RNA walk' a novel approach to study RNA-RNA interactions between a small RNA and its target.

In this study we describe a novel method to investigate the RNA-RNA interactions between a small RNA and its target that we termed 'RNA walk'. The method is based on UV-induced AMT cross-linking in vivo followed by affinity selection of the hybrid molecules and mapping the intermolecular adducts by RT-PCR or real-time PCR. Domains carrying the cross-linked adducts fail to efficiently amplify by PCR compared with non-cross-linked domains. This method was calibrated and used to study the interaction between a special tRNA-like molecule (sRNA-85) that is part of the trypanosome signal recognition particle (SRP) complex and the ribosome. Four contact sites between sRNA-85 and rRNA were identified by 'RNA walk' and were further fine-mapped by primer extension. Two of the contact sites are expected; one contact site mimics the interaction of the mammalian Alu domain of SRP with the ribosome and the other contact sites include a canonical tRNA interaction. The two other cross-linked sites could not be predicted. We propose that 'RNA walk, is a generic method to map target RNA small RNAs interactions in vivo.